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SPHERICAL  SHOCK  VAVF 
Richard  Latter 

The  RAJTD  Corporation,  8anta  Monica,  California 

ABSTRACT 

Tho  point -source ,  spherical  shock  wave  moving  into  a  constant  density, 
7-lav  gas  is  considered  in  the  limit  of  infinite  shock  strength  Tror.  tm 
point  of  view  of  the  Rlchtnryer-von  Neumann  viscosity  technique.  A  simi¬ 
larity  solution  of  this  problem  is  shown  to  exist  anl  is  obtained  for 
various  boundary  condition#  with  y  ~  l.k,  Tha  solutions  are  obtained  ana¬ 
lytically  in  that  part  of  the  flov  field  not  involving  viscosity,  and  numer¬ 
ically  in  the  other  porta  of  the  flov  field.  It  in  found  that  where a*  all 
discontinuities  of  the  physical  parameters  are  removed  by  the  viscosity  there 
remain  discontinuities  in  the  slope*  of  these  parameters  at  the  shock -front. 

It  is  indicated,  moreover,  that  the  complete  flov  field  depends  upon  tha  form 
and  magnitude  of  the  viscosity. 

I.  INTRODUCTION 

Tha  existence  of  shock  waves  in  n  hyd  dynamic  flov  field  introduces 
fraa-boundary  discontinuities  into  the  physical  parameters  of  the  syste;:.. 

Such  discontinuities  cause  considerable  analytic  as  well  as  numerical  compli¬ 
cations  for  the  trautment  of  hydrodynamic  problta*.  A  means  for  avoiding 
these  difficulties,  particularly  for  numerical  calculations,  hAS  been  developed 
by  Rlchtwyer  and  von  Neumann^.  They  observed  that  the  addition  of  a 


(l)  R.  D.  Richtcyer  and  J.  von  Neumann,  Journ.  of  Appl.  Phys.  21  2jf  ( 19l»0). 


particular  viscosit.. -like  tern  into  the  hyurodyuo.  ic  equations  could  lead  to 
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continuous  shock-flow  field*  in  which  the  discontinuities*  * t  the  shock  wavs 
were  removed  and  replaced  by  a  region  in  which  the  physical  parameters 
changed  rapidly,  but  smoothly.  Moreover,  the  change  of  the  physical  param¬ 
eters  through  this  transition  region  was,  for  piano  flow  fields,  in  sgres- 
aent  with  the  Rank ine -Eugenio t  relations. 

Riehtoysr  and  von  Keuacnn  obtained  an  analytic  solution  for  flow  field 
Including  the  artificial  viscosity  only  for  a  plane  shock  wave.  In  the  present 
discussion,  it  is  shown  (using  the  viscosity  technique)  that  in  the  limit  of 
infinite  shock  strength,  a  point-source  spherical  shock  flow  has  s  similarity 
solution  for  a  7 -lav  gas.  The  latter  has  been  obtained  by  numerical  and 
analytic  integration  of  the  ordinary  differential  equations  resulting  from  the 
similarity  of  the  flow. 

An  important  consequence  of  the  similarity  solution  described  here  is 
that  it  reveals  tnat  though  the  discontinuities  of  the  physical  parameters  sre 
removed  from  the  flow  field,  the  rates  of  charge  of  these  parameters  have  dis¬ 
continuities.  Moreover,  the  flow  field  itself  is  modified  by  the  viscosity, 

(2) 

approaching  the  Taylor'  atrong-shock  field  only  in  the  limit  of  vanishing 

(2)  G.  I.  Taylor,  Proc.  Roy.  Soc.  201,  159  (1950). 

viscosity. 

ii.  yoRWUTioa 

When  the  viscosity  teric  suggested  by  Ricbtmysr  and  von  tteumann  is  in¬ 
cluded  into  th#  hydrodynarlc  equations  for  s  spherically  symmetric  flow  in  a 
7-lav  medium,  tise  hydrodynamic  equations  in  Bulerian  form  becone 
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where  p  Is  the  static  pressure ,  p  the  density  of  the  aediura,  u  the  material 
velocity,  q  the  artificial  viscosity,  and  (r,t)  the  space-tine  Euler ian  coor¬ 
dinates.  The  caloric  equation  of  state  of  the  medium  is  assumed  to  be 

«  »  p/(r-l)p  ,  (4) 


where  <  is  the  Internal  energy  per  unit  mass  and  y  >  1.  The  adiabatic  equa¬ 

tion  of  state  is  assumed  to  be 


p/p7  -  o(8) 


(5) 


where  o(*)  is  a  function  of  the  entropy  8  only.  The  quantity  q  in  Kqs.  (l)  - 
(3)  is  the  viscosity  term.  That  dependence  of  q  on  p,p,  and  u  which  prohibits 
discontinuities  in  the  physical  parameters  is  not  uniquely  prescribed  by 
Richtrayer  and  von  Heumann.  The  fora  for  q  adopted  for  the  present  discussion 
is  consistent  with  their  requirements  and  Is 


q 


/ 

chi  ‘  !  chi 

Sr.  1$ 
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(6) 


The  arguments  presented  here  are  not  restricted  to  the  fora  of  q  in  JSq.  (6), 

but  can  be  generalized  to  include  other  possibilities.  However,  ao:ae  of  the 

detailed  results  depend  upon  the  specific  form  of  q  which  is  used}  in  particular 

not  all  forms  for  q  admit  of  similarity  solutions. 

Tha  possibility  of  a  similarity  solution  of  Eqs.  (l)  -  (3)  and  (6)  is 

(2) 

suggested  by  the  calculations  of  Taylor  on  the  point -source,  strong -shock 


f 


relation  to  the  shock  position,  pQ  is  an  arbitrarily  prescribed  constant  vith 
the  dimensions  of  density,  and  a  and  3  art  constants  yet  to  be  determined.  To 


simplify  the  subsequent  calculations  the  result  3  «  2a,  which  nay  be  easily 
derived  from  the  arguments  to  follow,  is  assumed  without  proof.  The  expres¬ 
sions  in  Sqa.  (7)  -  (9)  are  now  substituted  into  S^s.  (l)  -  (3)  and  (6), 
using  the  relationships,  which  result  from  the  change  of  independent  variables 
from  (r,t)  to  (x,t), 

e 

3  d  xR  i 

3t  -  3t  -  r  Sc  ' 

5  1  b 

Sr  ‘  T  St  ' 

where  ft  -  dfl  dt.  Equation  (l)  then  beco-'s  in  terms  of  the  new  variables 


JpR  -  +  #  +  LJJL. 
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where  the  prize  Indicates  differentiation  with  respect  to  x.  Tbs  quantity  g 
is  related  to  q  by  the  equation 


p~  "  ♦  *  i  L\  - 1  1=*  *4^  > 

po  2R  V 


,2a 


(U) 


\ 


which  or isos  directly  from  Bq.  (6)  using  Iqs.  (7)  -  (9).  If,  now,  it  is 
assumed  that 


A,  a  constant  , 


(12) 


then  Eq.  (10)  becomes  an  ordinary  differential  equation  in  the  independent 
variable  x,  namely, 

-  A(o^  +  xf$* )  -f  "*  *  0  .  (13) 

Quite  similarly,  Bqs.  (2)  and  (3)  become,  respectively, 

-Axt*  ♦fSf*  +  t*S*+2  &  •  0 


and 


-  A  2a  f  +  (^-Ax)f  -  (p-Ax)t'  »  0  . 


(14) 


(15) 


The  constant  a  is  determined  by  requiring  a  constant  energy  in  the  hydrodynamic 
flov.  If  E  Is  the  total  energy,  then 


00 


4*  * 

w 


r^dr  (pu^/2  ♦  pT'-l)  , 


ao 


which  becomes  in  terms  of  the  present,  similarity  variables 


00 


8  ■  ^  js  %  [ 


x*dx  iwS  2  +  f,fr-l)  )  . 


(16) 


Since  R  is  a  function  of  time,  this  equation  requires  that  a  »  3,  2. 


If  nev  diaeniionAe**  quantities 
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f  a  f 

F  *  —x  ,  G  *  "s-  f  and 

A  A 


1 


are  introduced  into  Sqt.  (11)  and  (13)  -  (l6) ,  these  equations  read 


o.£  -  f )  , 


/  3  •  •  '  -  •  F'+0«  . 

+—  ' 

I 

(f  -  x)t  1  +  U'  +  2  t  ~  -  0  , 

i  I 

-  3F  ♦  (:  -  x)F*  -  (  <5  .  x)  ♦ 


*  -  0  , 


and 


ao 

X  -  1*  p0  A2  f  x2  dx(t  |  2/2  ♦  F/7-1)  . 


(11* ) 

(13*) 

(1M) 

(15*) 

(16*) 


The  existence  of  a  solution  of  Equations  (13*)  -  (15*)  for  a  diverging  flow 
from  a  point  source  with  a  shock -discontinuity  at  x  =>  1  has  been  demonstrated 
by  Taylor.  But  for  such  a  diverging  flow  ^  >0,  so  that  G  «*  0,  and  therefore 

the  Taylor  solution  constitutes  also  a  solution  of  the  present  problem  in  the 
region  from  x  «  0  to  x  »  1,  where  R  is  identified  with  the  position  of  the 
Taylor  shock -front.  In  the  Taylor  solution  the  quantities  F,  |  and  t  at 
x  *  1  are  consistent  with  the  Rank ine -Hugoniot  conditions  for  a  strong  flfaock. 

The  presence  of  viscosity  removes  the  discontinuities  in  the  physical  quanti¬ 
ties  and  consequently,  the  Rankine -Hugoniot  condition*  have  no  special  signifi¬ 
cance  ir.  the  viscosity  formalism.  However,  if  the  viscosity  is  to  leave  the 
flov  field  unaltered,  then  the  Rank ine -Bugoniot  conditions  should  apply  across 
the  shock -transit ion  region  wherein  the  viscosity  Is  effective.  In  a  similarity 


problem  of  the  present  type  it  is  clear  that  this  situation  la  not  possible 
sines  the  vidth  of  tbs  transition  region  vili  be  a  constant  with  respect  to 
the  similarity  variable  x  -  jp  .  The  transition  region  will  there  fort'  con¬ 
tain  an  amount  of  mass  given  by 

n(r,t)  »  4a  r2dr  p(r,t)  , 

which  in  terns  of  the  similarity  variables  is 
*(r,t)  *  4*pQR^  j  *  ♦  (*)  • 

The  limits  of  integration  are  fixed  by  the  width  of  the  transition  region 
and  depend,  consequently,  only  upon  x.  The  time  dependence  of  R  given  in 
Sq.  (12)  implies  thee  that  the  mass,  m{r,t),  within  the  transition  region 
varies  with  tine.  For  the  validity  of  the  Ranklne -Hugoniot  conditions  across 
the  transition  region,  m(r,t)  would  have  to  be  a  constant,  but  this  condition 
has  been  proved  to  be  violated.  This  impossibility  of  meeting  the  Rank in* - 
Hugoniot  conditions,  la,  of  course,  a  property  of  the  assumed  type  of  simi¬ 
larity  expressed  by  lqs.  (7)  -  (9)  as  veil  aa  of  the  specific  form  for  q  of 
Bq.  (6).  The  restrictions  on  the  form  of  q  leading  to  the  Rank in* -Hugoniot 
conditions  are  etill  not  known  in  general  and  certainly  need  further  explora¬ 
tion. 

Before  studying  the  behavior  of  the  differential  equations,  it  ig  con¬ 
venient  to  formulate  the  boundary  conditions  which  tliese  equations  oust  sat¬ 
isfy.  Hamely,  all  physical  quantities  must  be  bounded  and  continuous.  While 
the  artificial  viscosity  is  not  a  physical  quantity,  spirit  of  tiie 
Richtmyor-von  Neumann  formalism  requires  treating  It  as  such  and  therefore  q 
is  also  assumed  continuous.  For  the  presently  considered  diverging  flow- 
field,  the  physical  quantities  must  satisfy  the  additional  con.lir.ions  that  In 

•*  0,  j>  •  0,  F  >  0,  G  *>  0,  and  t  >  0,  and  iu  front  of  the 


the  neighborhood  of  x 
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flov  the  medium  Is  undisturbed  so  that  q  *  0,  7  >,  0,  C  »  0,  and  f  is  a  con¬ 
stant.  From  Iq.  (7)  it  is  readily  observed  that  the  similarity  of  the  flov 
requires  1  *  0  for  large  x  in  the  undisturbed  medium. 

The  treatment  of  Equations  (11*)  and  (13*)  -  (15* )  vill  be  divided  into 
the  consideration  separately  of  the  regions  where in  the  viscosity  is  absent 
and  wherein  it  ie  preeent.  The  region  without  viscosity  will  be  treated  first 
and  vill  be  solved  analytically.  The  region  with  viscosity  will  then  be  ex¬ 
plored  end  joined  on  by  numerical  integration  to  the  rest  of  the  flov -fie id. 


in,  th  ruN-rmj)  vitboot  viscosity 

from  the  boundary  conditions  for  the  diverging  flow-field,  it  foil  ova 
that  in  a  neighborhood  of  x  -  0,  the  gradient  of  the  velocity  la  positive  and 
therefore  the  viscosity  -conn  0  is  sero.  The  equations  defining  the  flow  then 
simplify  to 


« 

0  ♦  X  0 


'  * 

+  00 


♦  0 


(IT) 


e  w 

(o-x)f  *  ^  0  ♦  2  t ' 


(18) 


-  3F  ♦  (o  -  x)y’  -  2L  (9  -  x)t*  -  o  .  (19) 

These  equations  admit  of  a  straightforward  analytic  integration.  Thus  by 
dividing  Bq,  (l8)  by  * (o  -  x)  end  by  adding  the  result  to  Eq.  (19),  after 
dividing  the  latter  by  F(o  -  x), 


V 

f 


(y  -  l)  -£-  + 

1 


0 


This  equation  integrates  directly  to 

xc(x  -  0)  «  a  .  (20) 

A  second  integral  of  Iqe.  (17)  -  (19)  °ay  be  obtained  by  nwltiplying  Bq.  (17) 


vsr-i 
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with  and  by  adding  to  Bq.  (19)  multiplied  by  .  Simplifying  the 

result  vith  the  aid  of  Iq.  (l8)  gives 

r 

^  jr  p  •  (x~p)  ♦  •£-  ♦  o2)}  -  o 

Integration  of  thie  equation,  using  the  condition  0  *  o  at  x  *  0  yield* 

*  m  10 

T  P  -  (x  -  0)  t  0  )  ■  0 


(21) 


Alternatively,  thi*  equation  is  as  lanediate  consequence  of  the  conservation 
of  energy  if  this  conservation  is  expressed  in  integral  fori?..  Finally  the 
third  integral  Is  obtained  by  diract  substitution  of  8qs*  (20)  and  (21)  into 
Bqs.  (17)  -  (19).  Thus,  first  eliainate  F*  andf  9  from  Bq.  (17)  using  Bq*. 
(18)  and  (19))  then 


9 

1 


0  (x-o)  If  +  (3-27  j) 
7  -  (x  -  i)2  | 


(22) 


f 

Expressing  t  and  F  in  terras  of  p  and  x  gives  for  0 

3(27-1)  ?  + 


T  • 

P 


p  27(7*1) -  3(27-1)  *  ♦  3 
7(7+1)  (Jf-  2(7+1)  \  +  2 


X 


(23) 


This  is  a  standard  homogeneous  equation  and  readily  integrated  giving 


x  »  b 


?  1 

«p»  « 


(7-1)/ (27+1) 


(24) 


f. 

*  * 


(l3r2-77+12)/5(27+l) (37-1) 


'3&r 


Xqs.  (20), (21)  and  (24)  constitute  a  cccqOete  analytic  solution  of  the 
differential  equations,  Xqs.  (l?)  -  (19)*  Since  F  and  *  are  intrinsically 
positive,  the  constant  a  of  Bq.  (20)  it  positive. 
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The  behavior  of  the  flow-field  in  the  present  region  it  defined  by 
fixing  the  two  arbitrary  integration  constant*  a  and  b.  These  are  in  turn 
specified  by  requiring  definite  values  for  £  and  F  at  a  prescribed  position 
x.  It  is  clear  therefore  that  the  Iqs.  (17 )  -  (19)  describe  a  two  parameter 
farily  of  flow-fields.  Ho  further  restriction  of  these  parameters  is  possi¬ 
ble  before  discussing  the  region  containing  viscosity. 

Only  the  properties  of  the  present  solution  relevant  to  the  viscosity 
fornaliem  will  be  considered.  First ,  sines  |  «  0  at  x  «  0,  Bq.  (24)  requires 
that  i  *  i  at  x  -  0.  8olving  Bq.  (21)  for  J  gives 


F  7zl  2\ xrl  "  x  * 

♦  *  aT  ; 

1  1 


X  7 


and,  moreover,  froir  the  poeltlvlty  of  F  andt  ,  this  equation  above  that 

1  <  '  <  , 

—  m  «•  as  i  * 

7  x 

Combining  *q».  (24)  and  (25)  leads  to  the  restriction 


7  *  7  *  Mln(1’  5^r >  • 


Using  this  equstion  with  Xq.  (23)  it  is  shown  immediately  that 


6  >  0  ,  (27i 

for  all  x  ^  0  provided  that  the  physical  require  sent  7  >_1  is  met.  This 

?  1  1*1 

result  is  established  by  observing  first  that  at  x  •  0,  j  «  -  sad  p  »  y  >  0. 

If  then  j>*  is  to  violate  Bq.  (27),  there  must  exist  s  point  at  which  ^  ■  0. 
Burt  this  is  possible  only  if 


27(7-1) 


If- 


3(27-1)  -  ♦  3  -  0  . 


'‘4  *  '■  >4^’^ 
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Solving  this  equation  for  ~  yields 

t.  3(ar-l|^  . 

It  may  os  demonetrated  analytically  from  this  txpreasion  or  simply  by  plot  - 
t  lng  the  two  roots  as  a  function  of 7  that  far 7  >  0 

I 

9 

But  this  result  implies  that  for  no  Talus  of  -  in  the  allowed  rang®  sx- 

T  f 

pressed  by  If.  (26)  does  ^  -  0. 

,  %  T« 

Since  the  viscosity  tern  of  lq.  (11')  remains  aero  so  long  as  0  >  0, 
there  Is  no  way  far  the  viscosity  to  aster  except  by  the  Introduction  of 
discontinuities  in  the  slopes  of  the  physical  quantities.  This  indicates  a 
deficiency  of  the  vlseoelty  formalism  which  in  numerical  applications  assumes 
continuity  not  only  of  the  physical  quantities  theme# lvee  but  also  of  their 
elopes. 

IV.  THE  flOW-FULD  WITH  VISCOSITY 

From  the  above  arguments  it  has  been  established  that  Kqa.  (13’ )  *  (15') 
ire  solved  by  a  Taylor-type  solution  in  the  region  x  ^  0  wherein  0*0.  But 
la  order  to  Introduce  the  vlseoelty  it  has  been  proved  necessary  to  adc.it  dis¬ 
continuous  slopes.  The  position  at  which  the  discontinuities  ere  located  may 
be  prescribed  arbitrarily  as  x  *  1  since  this  affects  only  the  magnitude  of  R 
which  Is  at  present  arbitrary  and  fixed  only  after  specifying  the  constant  A. 

The  problem  remaining  Is  then  the  solution  for  the  region  x  ^  1  In  which  the 
vlseoelty  le  effective.  This  region  comprises  the  transition  from  the  contin¬ 
uous  flow-field  to  the  undisturbed  medium  In  front  of  the  flow.  For  this 
region  •  0  and  therefore  Iq.  (11')  may  be  written 


■  yi- 
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Nf 


(n»») 


At  ths  connection  point,  x  -  1,  G  -  0,  a*  required  by  continuity  of  the 

- 1 

physical  quantities  and  therefore  f  >  0  at  this  point.  This  in  turn  deter- 

• » 

mines  the  Magnitude  of  the  Juqp  in  the  slope  b  ,  which  combined  with  the 
differential  equations  prescribes  the  Juqps  of  the  slopes  of  7,  ♦  and  G. 

The  solution  of  Sqs.  (11**)  and  (13*)  -  (15*)  for  the  transition  region 
must  aeet  the  boundary  conditions  that  for  sufficiently  lsrge  x,  the  nedivaa 
is  undisturbed  so  that  b  *  0,  F  *  0,  G  »  0  and  t  is  a  constant,  and  for 
x  »  1  the  solution  agrees  with  the  solution  for  the  region  x  *  0  to  x  »  1, 
wherein  0*0.  It  nay  be  shown  quite  generally  by  a  eerles  expansion  for 
large  x  that  all  solutions  of  Bqs.  (11**)  and  (13*)  *  (15*)  which  extend  to 
infinity  be hare  asymptotically  as 


e 


and  t_f  Q 


(26) 


t 


where  ,  fQ  ,  gQ  ,♦  ,  andt  1  are  constants.  This  result  may  be  verified 

by  direct  substitution  into  the  basic  equations.  This,  of  course,  establishes 


that  the  medium  In  front  of  the  shock -transit ion  region  is  undisturbed  and, 
moreover,  is  at  aero  pressure,  which  laplies  aero  temperature .  However,  not 
every  solution  extends  to  infinity  nor  is  physically  admissible  since  4>  ,  F,  G 
and*  must  be  everywhere  positive  in  order  to  correspond  to  a  Aivorging  flow- 
field  with  positive  pressures  and  densities.  It  may  be  shown  that  the  physically 
admissible  solutions  would  extend  to  infinite  x.  for  admissible  solutions, 
the  values  of  the  quantities  $,7,0  and*  at  x  -  1  must  be  restricted.  To 
examine  this  restriction,  it  is  observed  first  that  certain  restrictions  on 
ti©s»  quantities  already  exist.  Hasaly,  0  -  0  at  x  ■  1.  Also 


(M  fl)  |  •  (i  •  |)  (~*j*  ♦  ^  )  » 


0 


t 


(29) 


or  at  x  -  1 
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T  $  -  (l  -  0)  {^j*  ♦  j  fo2)  -  0  .  (30) 

Sq.  (29)  follow*  from  Xqs.  (13*)  -  (15® )  and.  Iq.  (20)  by  a  derivation  com¬ 
pletely  samlogous  to  that  leading  to  Bq.  (21) }  and  Eq.  (30)  1*  a  coneequence 
of  Q  •  0  at  x  •  1.  finally,  it  1*  observed  from  Eqs.  (11M)  and  (13*)  *  (15*) 
that  all  solution*  with  tbs  sens  values  of  Y,  ♦  ,  0/ 1  and  0  at  1  »  1  lead  to 

th.  «r.  function.  £|*|  (  ^  ^  0(l).  slace  G  -  0  at  x  ■  X  and 

Sq.  (30)  relates  f/*  to  0  at  s  »  1,  th*  only  independent  degree  of  free don 
remising  for  the  e  ho  lee  of  the  quest  it  let  at  x  -  1  with  which  to  obtain  an 
admissible  solution  Is  the  specification  of  ©  at  x  *  1.  That,  in  fact,  a 
choice  of  o  cas  be  made  which  lead*  to  as  admissible  solution  has  not  been 
established  analytically.  However,  a  convincing  demonstration  of  the  exis¬ 
tence  of  such  a  eboice  is  provided  by  the  explicit  numerical  solution.*  die- 
ouesed  is  the  next  section.  It  aay  be  considered  proved,  therefore,  that 
®qs.  (H**)  and  (13*)  -  (15*)  admit  of  solution*  which  consist  of  a  contin¬ 
uous  hydrodynamic  flow  from  x  »  0  to  x  *  1,  and  which  exhibit  a  continuous 
transition  region  (except  far  slope  discontinuities)  from  the  flow  conditions 
at  x  »  1  to  an  undisturbed  medium  at  xero  pressure  for  large  x. 

V.  HMRJLT8  AID  PIflCUMIO X 

To  specify  the  diverging  flow  field  the  constants  a  and  b  of  Eqs.  (20) 
and  (2*0  (which  nr*  equivalent  to  o  (1)  and  F(l)/t(l))  and  the  constant  A  of 
Bq.  (12)  must  be  fixed.  The  energy  E  of  Eq.  (16#)  fixes  one  relationship 
amongst  these  constants  and  the  specification  of  the  density  of  the  undis¬ 
turbed  medium  fixes  another  relationship.  Finally,  the  requirement  of  a 
physically  admissible  solution  in  the  transition  region  fixes  the  last  constant. 
In  th©  limit  of  vanishing  viecosity,  1*0,  the  solutions,  of  course,  becone 
identical  with  the  familiar  strcag-ehock  flow  fields  with  a  Jump  discontinuity 


and  for  I  »  0. >9  it  is  1.137.  Eta  ttsa  of  the  viscosity  foraaliaa  is 
•ten  fro®  these  reties,  M  will  as  from  the  complete  solutions  in  Tables  I  - 
17,  to  introduce  considsr&bls  modification  of  the  flov-field  for  large  K,  but 
negligible  changes  for  saall  X. 

The  important  eoaeequsKie  of  the  present  discussion  is  the  quantitative 
estimate  of  the  influence  of  the  viscosity  on  the  flow  field.  Perhaps  as 
significant  is  the  proof  that  the  viscosity  technique  may  possess  the  defi¬ 
ciency  in  spherical  flow  fields  that  the  physical  parameters  have  discontin¬ 
uous  slopes  at  the  shock-transition  region,  even  though  these  parameters  are 
cocrtinuaus.  Thus  the  usual  numerical  Integration  of  Bqs.  (l)  -  (3)  and  (6), 
carried  out  by  replacing  these  equations  with  difference  equations,  cannot 
lead  to  a  solution  of  tbs  differential  equation  if  the  slops  diecontinuitlss 
are  disregarded.  It  would  appear  therefore  that  these  discontinuities  should 
lead  to  inaccuracies  in  rxaoarlcal  applications  of  the  viscosity  technique, 
which  ignore  them.  However,  free  3.  Brode* a ^  results  on  the  numerical 


(3)  H.  Brode,  to  be  published. 

integration  of  the  hydrodynamic  Sqs.  (1)  -  (3)  using  the  viscosity  technique 
for  a  spherical  explosion,  it  appears  that  any  error#  from  this  source  are 
probably  unimportant.  This  conclusion  should  be  qualified  somewhat  since  In 
Brode* s  integration,  the  radius  to  the  transition  region  was  large  and  the 
neighborhood  of  the  shock -front  was  approximately  planar  sc  that  any  discon¬ 
tinuities  would  be  small.  Moreover,  the  form  of  viscosity  used  by  Brode 
differed  from  the  present  one,  although  the  form  which  he  used  also  suffered 
fro®  act  leading  to  a  constant  m&a,  within  the  transition  region. 
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The  author  wishes  to  express  his  appreciation  to  Mr.  Jack  Little  and 
Miss  Ruth  Merrill  for  numerical  assistance. 
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Captions  for  Figuras 


Figure  1 
Figure  2 
Figure  3 
Figure  h 


Solution  for  K  «  0 
Solution  for  £  *  0.003^9 
Solution  for  X  «  0.03^9 
Solution  for  X  ■  0.3^9 
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